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ABSTRACT 

We analyse the evolution of the properties of the low-redshift Intergalactic Medium 
(IGM) using high-resolution hydrodynamic simulations that include a detailed chem- 
ical evolution model. We focus on the effects that two different forms of energy feed- 
back, strong galactic winds driven by supernova explosion and Active Galactic Nuclei 
(AGN) powered by gas accretion onto super-massive black holes (BHs), have on the 
thermo- and chemo-dynamical properties of of the low redshift IGM. We find that 
feedback associated to winds (W) and BHs leave distinct signatures in both the chem- 
ical and thermal history of the baryons, especially at redshift z < 3. BH feedback 
produces an amount of gas with temperature in the range 10 5 — 10 7 K, the Warm Hot 
Intergalactic Medium (WHIM), larger than that produced by the wind feedback. At 
z = the fraction of baryons in the WHIM is about 50 per cent in the runs with 
BH feedback and about 40 per cent in the runs with wind feedback. The amount of 
warm baryons (10 4 < T < 10 5 K) is instead at about the same level, ~ 30 per cent, 
in the runs with BH and wind feedback. Also, BH feedback provides a stronger and 
more pristine enrichment of the WHIM. We find that the metal-mass weighted age of 
WHIM enrichment at z = is on average a factor ~ 1.5 smaller in the BH run than for 
the corresponding runs with galactic winds. We present results for the enrichment in 
terms of mass and metallicity distributions for the WHIM phase, both as a function of 
density and temperature. Finally, we compute the evolution of the relative abundances 
between different heavy elements, namely Oxygen, Carbon and Iron. While both C/O 
and O /Fe evolve differently at high redshifts for different feedback models, their values 
are similar at z = 0. We also find that changing the stellar initial mass function has 
a smaller effect on the evolution of the above relative abundances than changing the 
feedback model. The sensitivity of WHIM properties on the implemented feedback 
scheme could be important both for discriminating between different feedback physics 
and for detecting the WHIM with future far-UV and X-ray telescopes. 

Key words: Cosmology: theory - Methods: Numerical - Galaxies: Intergalactic 
Medium 



1 INTRODUCTION 

Solving the problem of the missing baryons at low red- 
shift is one of the important goals of observational cos- 
mology (e.g.. iPersic fe Saluccil Il992l ; iFukugita et"afl 1 19981 : 
Fukugita fc Peebles! 120041 ). Since the pioneering work by 
Cen fe Ostrikerl (119990 ■ cosmological hydrodynamic simula- 
tion s provided a fund amental contribution in this field (see 
also lDave et aLlfeOQll ). Although the baryon budget is closed 
by observations at high redshift, z£, 2, simulations indicate 
that about half of the baryons at low-z should lie in a tenu- 



ous quite elusive phase, the so-called Warm-Hot Integalactic 
Medium (WHIM). The WHIM is predicted to be made of 
low density (nn ~ 10 -6 — 10~ 4 cm -3 ) and relatively high 
temperature (T ~ 10 5 — 10 7 K) plasma, whose main con- 
stituents are ionised Hydrogen and Helium, with traces of 
heavier elements in high ionization states. Simulations have 
also shown that gas at these densities traces the filamentary 
structures which define the skeleton of the large-scale cosmic 
web, and is heated by shocks from supersonic gravitational 
accretion onto the forming potential wells. 

Due to its low density, the collisionally ionised WHIM 
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should be characterised in emission by a very low surface 
brightness in the UV and soft X-ray bands, so that its 
detection lies beyond the capability of available detectors 
and shoul d await for instrumen tation of the next genera- 
tion (e.g.. lYbshikawa et alj 12004 ) ■ Claims for the detection 
of the WHIM t hrough soft X-ray emission have already bee n 
presented (e.g.. lZappacosta et al.|[2005l ; IWerner et al.l l2008). 
However, these detections are limited to special places, such 
as external regions of massive galaxy clusters or large fil- 
aments connecting such clusters, where the gas density 
reaches values higher than expected for the bulk of the 
WHIM. A more promising approach to reveal the pres- 
ence of the WHIM lies in the detection of absorption fea- 
tures in the spectra of background sources, at far ultraviolet 
(FUV) and soft X-ray energies, associated to atomic transi- 
tions from highl y ionized elements, using GRB (Gamma Ray 
Burst) spectra (jFiore et alj|2000l ; iBranchini et al.ll2009f ). or 
in the detection of Ly-a absorption from the tiny fraction 
of n eutral hydrogen rev ealed through FUV spectroscopy 
(see iRichter et all [2008, for a review). Although detec- 
tion of gas with T < 10 5 K has been obtained at FUV 
frequencies along a fairly large number of sightlines (e.g., 
iTripp et alj|2008l . and references therein), the situation is 
less clear in the soft X-ray band. At these energies, claims of 
WHIM detection in absorption have be e n presented by di f- 
ferent authors (e.g. iNicastro et al.ll2005l ; iBuote et al.ir2 009). 
although they are either con troversial (e.g. iKaastra et all 
2006; iRasmussen et all 120071 ) or have moderate statistical 
significance. 

Owing to the above observational difficulties in detect- 
ing and characterizing the physical properties of the WHIM, 
numerical simulations have played over the last years the 
twofold ro le to forecast its detectability, both in emis- 
sion (e.g., iRoncarelli et al. 20061; Ursino fc Galeazzi 20061) 



and i n absorption (e . g.. ICen et al. 2001 ; Kravtsov et alj 
120021 ; IViel et all [200I ; IChen et al.l 120031 ; IViel et all I200IE 

and to predict its observational properties, also as a func- 
tion of the physical processes included (see iBertone et al.l 
2008, for a review). For instance, feedback effects from 
star formation and accretion onto super-massive black holes 
(SMBHs) are expected to determine at the same time 
the process of galaxy formation and the physical prop- 
erties of the diffuse cosmic baryons. However, as of to- 
day, the mechanisms powering feedback are still poorly un- 
derstood and thus very difficult to model in a fully self- 
consistent way in cosmological hydrodynamical simulations. 
On the other hand, numeric al simulations from different 
groups (e.g.. |Cen fc Ostrikerl 



2006; Dave & Oppcnhcimcr 
lOppenheimer fc Pavel 1200*1 
Tescari et al.l 120091 ; IWiersma et al.l 120091 ) confirm that the 
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thermo- and chemo-dynamical properties of the WHIM and, 
more generally, of the Intergalactic Medium (IGM) at dif- 
ferent redshifts are sensitive to the nature and timing of 
the feedback. Indeed, metals are observed in the low den- 
sity IGM out to high redshift and it is likely that galactic 
winds, as observed in galaxies in the low redshift universe, 
are responsible for the IGM metal enrichment. Furthermore, 
SMBHs are also expected to have an impact on the gas dis- 
tribution in galaxies and to blow ejecta out to large dis- 
tances, especially during galaxy mergers. 

In this paper, we will present an analysis of an extended 
set of cosmological hydrodynamical simulations with the 



purpose of investigating the effect that both galactic winds 
powered by supenova (SN) explosions and energy feedback 
from accretion onto SMBHs have in determining the ther- 
mal and chemical properties of the IGM. Although the main 
focus of the analysis is on the low-redshift IGM, we will also 
discuss how different feedback mechanisms leave their im- 
print on the evolutionary properties of diffuse baryons since 
z ~ 4. Our simulations are bas ed on the chem o-dynamical 
version of the GADGE T- 2 code (|Springell2005l) presented by 
iTornatore et al.l (|2007l ). which follows the production of dif- 
ferent chemical species by accounting for detailed yields from 
Type-la and Type-II SN (SN-Ia and SN-II hereafter), as well 
as from intermediate and low mass stars in the thermally- 
pulsating asymptotic giant branch (TP-AGB) phase, while 
also accounting for the mass dependent life-times with 
which different stellar populations release their nucleosyn- 
thetic prod ucts. The model of galactic winds is that in- 
troduced bv lSpringel fc Hernquistl (|2003l ). that we consider 
both in its original version and in a version based on assum- 
ing that winds are never hydrodyn amically decoupled from 
the s urrounding medium (see also iDalla Vecchia fc Schavel 
120081 ). As for the B H feedback, we adopt t h e models origi- 
nally introduced by |Pi Matteo et al.l |2005l) ; ISpringel et alj 
(2005). Besides investigating the effect of different feedback 
mechanisms, we will also consider the impact of changing 
the stellar initial mass function (IMF) on the resulting en- 
richment pattern of the IGM. Although the results of this 
paper have important implications for the detectability of 
the WHIM in view of next-generation X-ray missions, we 
defer an observationally oriented analysis of our simulations 
to a future work, aimed at investigating in detail how dif- 
ferent instrumental capabilities will be able to characterize 
the WHIM properties at the level required to discriminate 
between different feedback models. 

The plan of the paper is as follows. In Section [2] we 
describe the hydrodynamical simulations used and the feed- 
back mechanisms adopted. Section [3] describes the main re- 
sults of our analysis in terms of global gas properties (evolu- 
tion of different phases over redshift), epoch of enrichment of 
gas particles, metallicity distributions at z = as a function 
of overdensity and temperatures, and evolution of different 
chemical elements. We provide a final discussion of our re- 
sults and draw our main conclusions in Section [4] In the 
following, we will assume the v alues of the solar metallicity 
as reported by lAsplund et al.l (120051 ). with Z/X = 0.0165 
(X: hydrogen mass; Z: mass contributed by all elements 
heavier than Helium). 



2 HYDRODYNAMICAL SIMULATIONS 

Our simulations were carried out using a version of t he paral- 
lel hy drodynamical TreePM-SPH code GADGET-2 (|Springell 
120051 ). which includes a detailed implementation of th e 
chemical enrichment as described in ITornatore et all l|2007l ). 
The initial conditions are generated at redshift z — 99 
in a cosmological volume with periodic boundary condi- 
tions filled with an initially equal number of dark mat- 
ter and gas particles. Th e matter power spectrum i s gen- 
erated using CMBFAST (|Seliak fc Zaldarriagal Il996l ) for a 
fiat ACDM model, with cosmological param eters consistent 
with the recent findings of WMAP year-5 (|Komatsu et al.l 
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l2009h : SJom = 0.24, Q a = 0.76, n oh = 0.0413, for the den- 
sity parameters contributed by total matter, cosmological 
constant and baryons, n s = 0.95 for the primordial spec- 
tral index, Hq = 73 km s _1 Mpc -1 and erg = 0.8 for the 
normalization of the power spectrum. Besides performing 
simulations within boxes of 37.5 ft -1 comoving Mpc using 
2 x 256 3 gas and dark matter (DM) particles, we also run a 
few simulations within a twice as large box of 75 ft _1 Mpc on 
a side, with 2 x 512 3 particles so as to keep resolution con- 
stant. Such larger boxes will allow us to keep under control 
any effect related to box size. The resulting mass of the gas 
particles is thus m gas ~ 3.6 x 10 7 h~ 1 MQ. In these runs, the 
gravitational softening is set to 6 = 7.5 ft -1 comoving kpc 
above z = 2, while at z < 2 it is set to e = 2.5 ft -1 phys- 
ical kpc. As for the B-spline softening length used for the 
computation of the SPH forces, the lowest allowed value was 
set to half of the gravitational softening. In order to address 
resolution effects, we also carried out one simulation in the 
smaller box using 2 x 400 3 particles. In this case the mass 
of the gas particles is m gas — 9.4 x 10 6 ft -1 Mq, with all the 



softening lengths rescaled according to m, 



-1/3. 



4.8 ft" 1 



comoving kpc at z > 2 and e = 1.6 ft 1 physical kpc at lower 
redshift. All these simulations are evolved to z = 0. 

We point out that both a large volume and high resolu- 
tion are necessary in order to correctly describe the physical 
properties of the WHIM. Indeed, shocks driven by the col- 
lapse of large-scale structures are important, while dense en- 
vironments able to trigger and power the feedback processes 
need to be resolved. Our simulations of different resolutions 
and/or box-sizes will be denoted by the tuple (linear size of 
the box in comoving ft _1 Mpc, number of DM particles 1 ^ 3 ). 
For instance, we will indicate with (75,512) a simulation 
within a box of 75ft _1 Mpc and containing 512 3 DM parti- 
cles (see Table 1). 

Radiative cooling and heating processes are followed for 
a primordial mix of hydrogen and helium and we include the 
dependence of cooling on me tallicity by adopting t h e coo l- 
ing rates from the tables by ISutherland fc Dopital (|l993h . 
We also include the effect of a spatially uniform redshift - 
dependent Ultra Vi olet Background (UVB) produced by 
quasars, as given bv lHaardt fc Madaul (| 1996T ) . with helium 
heating rates multiplied by a factor 3.3 in order to better fit 
observational constraints on the temperature evolution of 
the IGM. This background gives natur ally a hydrogen ion - 
ization rate F_i 2 ~ 0.8 at z = 2 - 4 (jBolton et al.ll2005h . 
which is in broad agreement with the observations. We adopt 
the effective model of star for mation from a multiphase inter - 
stellar medium introduced bv lSpringel fc Hernquistl <|2003h . 
In this model, each gas particle whose density exceeds a 
limiting threshold value is assumed to contain a hot and a 
cold phase, the latter providing the reservoir of star forma- 
tion. The two phases coexist in pressure equilibrium with 
their relative fractions being computed according to the lo- 
cal conditions of density and temperature. In our simulations 
we assume that the density threshold for a gas particle to 
become multiphase and, therefore, star forming, is tih = 0.1 
cm -3 in terms of the number density of hydrogen atoms. 
Star forming gas particles are then assumed to spawn col- 
lisionless gas particles a ccording to the st ochastic scheme 
originally introduced by iKatz et al.l (|l996h . We allow each 
gas particle to produce up to three generations of star parti- 



cles, each having therefore a typical mass of about one third 
of the initial mass of the gas particles. 

Our simulations have been carried out with the chemo- 
dynamical version of the GADGET-2 code described by 
iTornatore et al.l (|2007l 'l . The included model of chemical evo- 
lution allows us to follow the production of six different 
metal species (C, O, Mg, S, Si, Fe) from Type-II and Type 
la supernovae (SNII, SNIa), along with low and intermedi- 
ate mass stars (LIMS) in the thermally-pulsating asymptotic 
giant branch (TP-AGB) phase. Besides including different 
contributions from SNII, SNIa and LIMS, we also include 
the effect of the mass-dependent time delay with which 
different stellar populations release metals. Specifically, we 
adopt the lifetime function given by Padovani fc Matteuccl 
(|l993h . We us e the stellar yields byrThielemann et all (|2003l ) 
for SNIa, by IWooslev fc Weaver] dl995h for SNII and by 
Ivan den Hoek fc Groenewegenl |l997h for LIMS. The mass- 
range for the SNII is M > 8Mq , while for SNIa arising from 
binary systems in the mass range it is 0.8 < M/Mq < 8, 
with a binary fraction of 10 per cent. Finally, we use three 
distinct stellar initial mass functions (IMFs): a Salpeter 
(1955), a Kroupa (2001) and an Arimoto-Yoshii (1987) IMF. 
Our reference choice is the functional form proposed by 
Kroupa (2001) which adopts a multi-slope approximation, 
ifi{m) oc m~ y with y — 0.3 for stellar mass m < 0.5 Mq, 
y = 1.2 for 0.5 M Q < m < 1 Mq, and y = 1.7 for 
m ^ IMq. Both t he IM F by ISalpeterl (119551 ) and that by 
lArimoto fc Yoshiil (|l987h have instead a single slope, with 
y — 1.35 and y — 0.95, respectively. Our model of chem- 
ical evolution also includes stellar mass losses, which are 
self-consistently computed for a given IMF and life-time 
function. This means that a fraction of the mass of the 
star particles is restored as diffuse gas during the evolu- 
tion an d distributed to the su rrounding gas particles. We 
refer to ITornatore et all (|2007T ) for a detailed description of 
our implementat ion o f chemo-dynamics in GADGET-2, to 
ISaro et all l^OOrj ) and lFabian et al.l J2008I) for applications 
to simulations of galaxy clusters and to Tescari et al.l (120091 ) 
for a description of the global IGM properties around DLAs 
at z > 2. 

In the two following subsections we will summarize the 
main features of the two different feedback schemes adopted 
in this paper. 



2.1 Galactic outflows 

The implem entation of galactic o u tflows is extensively de- 
scribed by ISpringel fc Hernquistl (|2003h . In this model, 
winds are assumed to blow from star forming regions with 
a mass loading rate M w proportional to the star formation 
rate according to M w = r]M+. Star- forming gas particles 
are then stochastically selected to become part of a blowing 
wind, with a probability which is proportional to their star 
formation rate. In its original implementation, whenever a 
particle is uploaded to the wind, it is decoupled from hy- 
drodynamics until the density of the surrounding gas drops 
below a given limiting value. This allows the wind particle 
to travel 'freely' up to few kpc until it has left the dense star- 
forming phase, without directly affecting it. As a protection 
against decoupling a wind particle indefinitely in case it gets 
'stuck' in the ISM of a very massive galaxy, the maximum al- 
lowed time for a wind particle to stay hydrodynamically de- 
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coupled is set to idee — lvr/v w> where we fix Z w = 10/i _1 kpc 
in our reference case, while u w is the wind speed. As for the 
limiting density for hydrodynamic decoupling of gas parti- 
cles, it is set t o 0.5 in units of the th r eshol d for star forma- 
tion. Unlike in ISpringel fc HernquistJ (|2003h . we decide here 
to fix the velocity of the winds to the value v w = 500 km s _1 , 
instead of fixing the fraction of the SNII energy powering 
galactic ejecta. For the efficiency of the wind mass loading, 
we adopt r\ = 2. In sum, four parameters fully specify the 
wind model: the wind efficiency rj, the wind speed t> w , the 
wind free travel length Z w (whose detailed value is unimpor- 
tant) and the wind free travel density factor. 

In order to verify the effect of the hydrodynamic de- 
coupling of wind particles we also perform a simulation 
without such a decoupling, tdec = 0, i.e. keeping the 
particles always hydro dynamically coupled. We note that 
iNagamine et al.l (|2007h showed that global DLAs properties 
are r elatively insensitive to the valu e of i w . On the other 
hand. lDalla Vecchia fc Schavd (|2008h pointed out that keep- 
ing winds always hydrodynamically coupled has a significant 
effect on the evolution and star formation in simulations of 
isolated disk galaxies. 

We stress that our wind model is not the only possible 
wind implementation and others could be adopted such as 
those based on mo ment um-driven winds as su ggested by 
iMurrav et ail <|2005h and lDave fc Qppenheimer] (j2007l ). 



each neighbouring gas particle a probability of contribut- 
ing to the accretion, which is proportional to the SPH ker- 
nel weight computed at the particle position. In the scheme 
described above, this stochastic accretion is used only to 
increase the dynamic mass of the BHs, while their mass en- 
tering in the computation of the accretion rate is followed 
in a continuous way, by integrating the analytic expression 
for Mbh- Once the amount of energy to be thermalised is 
computed for each BH at a given time-step, this energy is 
then distributed to the surrounding gas particles using the 
SPH kernel weighting scheme. 



To sum up, we run simulations with the lKroupal (1200 ll ) 
IMF by turning off galactic winds (NW) , with galactic winds 
(W), with galactic winds always hydrodynamically coupled 
(CW) and with black hole feedback (BH). We summarize 
in Table 1 the simulations analysed in this paper. Overall, 
our simulation set allows us to address the effect of chang- 
ing: (a) box size (by comparing W37,256 with W 75, 512 and 
Way37,256 with Way 7 5,5i2), (b) IMF (by comparing W37.256 
with Way37,256 and with Ws37 i 25e), (c) resolution (by com- 
paring W37.256 with W37,40o) and (d) nature of the en- 
ergy feedback (by comparing W37,256 with NW37,256, with 
CW37,256 and with BH 3 7 i2 56). 



2.2 BH feedback 

We also include in our simulations the effect of feed- 
back energy from gas accretion onto super-massive black 
holes (BHs), followin g the s cheme originally introd uced by 
jSpringel et al.ll2005l . see also lDi Matteo et al.ll2005T ). In this 
model, BHs are represented by collisionless sink particles ini- 
tially seeded in just resolved DM haloes, which subsequently 
grow via gas accretion and through mergers with other BHs 
during close encounters. Every new dark matter halo, iden- 
tified by a run-time friends-of-friends algorithm, above the 
mass threshold Mth = 10 hT 1 Mq, is seeded with a central 
BH of initial mass 10 5 H~ 1 Mq, provided the halo does not 
contain any BH yet. Each BH can then grow by local gas 
accretion, with a rate given by 

Mbh = min (A/ B ,M E dd) . (1) 

Here Mb is the accretion rat e estimated with the Bondi- 
Hoyle-Lyttleton formula (e.g.. lBondilll952T ). while MEdd is 
the Eddington rate. The latter is inversely proportional 
to the radiative efficiency e r , which gives the radiated en- 
ergy in units of the energy assoc iated to the accrete d 
mass: e r = L t /(Mbhc 2 ). Following ISpringel et al. I (|2005h . 
we use e r = 0.1 as a reference value, which is typical for 
a radiatively efficient accr etion onto a Schwarzschild BH 
l|Shakura fc Sunvaevlll"973h . The model then assumes that 
a fraction et of the radiated energy is thermally coupled to 
the surrounding gas, so that -Bf ee d = e Y etM-e,Yt<? is the rate 
of the ene rgy released to h e at the surrounding gas. Using 
et ~ 0.05, iDi Matteo et al.l (|2005h were able to reproduce 
the observed Mbh — o relation between bulge velocity dis- 
persio n and mass of the hoste d BH (see also ISiiacki et al.l 
feoOSt h lpTMatteo et all (120081 1 ). Gas particle accretion onto 
the BH is implemented in a stochastic way, by assigning to 



3 RESULTS 

In Figure [T] we show the star formation rates (SFR) for the 
W37.256, NW37,256, CW3 7 ,256 and BH37 2 56 sim ulations along 
with observational results by [Hopkins (2004). At 3 the 
NW and BH runs behave very similarly. This is not sur- 
prising since BH feedback is not effective until a sufficiently 
large number of DM haloes, massive enough to host a seed 
BH, is numerically resolved. After that, gas accretion takes 
place in BHs with a subsequent release of thermal energy. 
Once BH feedback becomes efficient, star formation is sud- 
denly suppressed by the expulsion of hot gas. Kinetic feed- 
back by galactic winds is relatively more efficient at high 
redshifts in depriving relatively small haloes of the star- 
forming gas, while it fails in regulating the star formation 
within massive haloes at later epochs. We note that hydro- 
dynamically coupled winds (CW) tend to produce a slightly 
higher SFR around the peak of the star formation history. 
This is due to the fact that coupled winds tend to ther- 
malize their kinetic energy at smaller distances from star 
forming and dense regions. The cooling time of this gas is 
correspondingly shorter, which causes a larger fraction of the 
thermalized wind energy to be radiated away, thus reducing 
the feedback efficiency. Both runs including winds reproduce 
the observed behaviour of the SFR at high redshift, z£, 2, 
while they tend to produce too high star formation at low 
redshift. On the contrary, the run with BH feedback has 
a too high star formation at z > 3, while it recovers the 
observed SF level at z'S 2. This illustrates once again the 
different role played by feedback related to star formation 
and to BH accretion: while the former is efficient since early 
times in regulating gas cooling within small haloes, the latter 
sets in at relatively lower redshift to quench star formation 
within recently formed massive haloes. 
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Run 


Box size 


N l/3 

D M 


IMF 


Feedback 


W37,256 


37.5 


256 


Kroupa 


Winds; v w = 500 km s -1 


* * iO,Qlz 


75.0 


512 


Kroupa 


Winds: w w = 500kms~^ 


W37 : 400 


37.5 


400 


Kroupa 


Winds; v w = 500 km s -1 


Way37,256 


37.5 


256 


Arimoto-Yoshii 


Winds; v w = 500 km s -1 


Ws 37i 256 


37.5 


256 


Salpeter 


Winds; d w = 500kms _1 


NW37,256 


37.5 


256 


Kroupa 


No feedback 


CW37,256 


37.5 


256 


Kroupa 


Coupled winds; n w = 500kms _1 


61137^56 


37.5 


256 


Kroupa 


Black Hole feedb., no winds 



Table 1. Summary of the different runs. Column 1: run name; column 2: comoving box size (units of h Mpc); column 3: number of 
DM particles; column 4: stellar initial mass function (IMF, see text); column 5: feedback included (see text). 
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2 4 6 8 10 12 14 
redshift 

Figure 1. Star formation rate density as a function of redshift in 
units of M0yr —1 Mpc~ 3 for W37,256 (winds, continuous blue line), 
NW37 i 256 ( no feedback, dashed green line), CW37 : 256 (coupled 
winds, dot-dashed orange line) and BH 3 7 i 256 (black hole feed- 
back, triple dotted-dashed line) runs. Observationa l data points 
represented as full circles are taken from iHopkind l|2004r i. while 
the empty circles are from Bouwcns et al. 



3.1 Global gas properties 

We show in Figure [2] the density-temperature phase di- 
agr ams for the g as in the four simulations based on 
the iKroupal (|200lT ) IMF and including different feedback 
schemes. The four panels are for the simulations without 
feedback (NW 3 7,256, top left), with feedback associated to 
galactic winds (W37,256, top right), with galactic winds al- 
ways hydrodynamically coupled (CW37,256, bottom left) and 
with BH feedback (31437,256, bottom right). Each phase di- 
agram is color-coded according to the mass in metals which 
is associated to the gas particles belonging to each two- 
dimensional bin in the p gas -T plane, with brighter regions 
corresponding to a larger metal mass. As such, these plots 
convey information on the effect that different feedback 



schemes have on the way in which metals are distributed 
in density and temperature. Here and in the following we 
do not include in the baryon budget inventory the inter- 
stellar medium (ISM), i.e. for gas having density larger than 
the density threshold for the star formation (see Sect. 2). We 
note that only a tiny fraction of the gas lies in this phase 
at all redshifts. Furthermore, while our sub-grid model for 
star formation should provide an effective description of the 
ISM, it is not obvious that this description provides a fully 
realistic description of its temperature. For this reason, our 
gas mass and metal budget will not include the contribution 
from star-forming gas. 

The effect of feedback is quite strong for gas in the 
temperature range 10 5 -10 7 K (corre sponding to the WHIM 
phase that we will define below; e.g. ICen fc Ostriker1l2006l ). 
Strong galactic winds increase the amount of metals car- 
ried by gas with temperature in the range 10 5 — 10 6 K and 
overdensityQ <5b — 1-10, when compared to the no feedback 
case. This is a consequence of the fact that winds are loaded 
with gas particles which were star forming and, as such, have 
been heavily enriched. The effect is even more dramatic in 
the presence of BH feedback: even regions below the mean 
density are enriched with metals up to rather high temper- 
atures (~ 10 7 K) when compared to the corresponding W 
and CW simulations. This counter-intuitive result follows 
from the fact that BHs' release of large amount of energy 
in a relatively short time interval, around z ~ 3. This sud- 
den energy release turns out to be much more efficient than 
winds to heat metal enriched gas to high temperatures, thus 
displacing it from the haloes of star-forming regions to low- 
density regions. Once brought to high entropy by BH feed- 
back, this enriched gas is prevented from re-accreting into 
collapsed haloes at lower redshift. Gas in photo- ionization 
equilibrium at T ~ 10 4 K contains a comparable amount 
of metals in the three runs which include feedback. In the 
NW case the amount of metals present in the gas around 
the mean density is an order of magnitude lower than in the 
runs including winds. Indeed, in the NW simulation, most 
of the enriched gas remains at high density, instead of being 
transported away from star-forming regions. The short cool- 
ing time of this enriched gas causes its selective removal from 
the diffuse phase into the stellar phase. We note also that rel- 
atively cold gas in very dense regions, p gas — (10 3 — 10 4 )pb, 

1 Here and in the following we denote with 5b the overdensity 
of gas with respect to the mean cosmic baryon density p"b, i.e. 
s b = Pb/pb - 1. 
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Figure 2. The p gas — T phase diagrams at z = for simulations based on the iKroupal j200lf) IMF. by varying the feedback scheme, 
color coded according to the metal mass fraction (see vertical bar): NWsy^ss r u n with no feedback (top left), W37,256 run with galactic 
winds (top right), CW37,256 run with galactic winds always hydrodynamically coupled (bottom left) and BH37,256 run with BH feedback 
(bottom right). 



is more enriched with metals in the NW, CW and W simu- 
lations than in the BH one. This demonstrates the efficiency 
of BH feedback in displacing highly enriched gas outside the 
core regions of virialized haloes, while correspondingly in- 
creasing the enrichment level of gas at lower density. As for 
the comparison between coupled and decoupled winds (W 
and CW runs, respectively), the latter has more enriched gas 
at densities approaching the threshold for the onset of star 
formation. This is the consequence of the hydrodynamic cou- 



pling of winds which causes wind particles to remain more 
confined in the proximity of star-forming regions. 

From a qualitative inspection of Fig. [2] we draw the 
following conclusions: i) galactic winds have a large impact 
on the metal enrichment of the IGM at z = both for 
underdense-cold and mean density-hot regions of the T — 
Pgas plane; ii) AGN feedback is more efficient than winds 
in enriching the warm-hot gas at T = 10 5 — 10 7 K and 
relatively low overdensity, 8b 50; iii) correspondingly, very 
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dense regions, with overdensity 5b > 10 4 , are less enriched 
with metals in the BH simulations than in the other runs. 



3.2 Redshift evolution of the different phases 

In order to make our results directly comparabl e with 
those from previou s analyses (e.g., iDave et al.l l200ll ; 
ICen fc Ostrikeill2006f ). we define four different phases for 
the baryons in our simulations: a WHIM phase, for gas par- 
ticles with temperature in the range 10 5 — 10 7 K; a warm 
phase, for gas particles with T = 10 4 — 10 5 K; a hot phase, 
for gas particles with T > 10 7 K; a condensed phase, which 
includes all the baryonic mass associated with star particles. 

In Figure |3j we show the redshift evolution of the mass 
fraction in these four different phases, comparing the ef- 
fect of changing feedback scheme. The W37.256, NW37,256, 
CW37,256 and BH37,256 runs are represented by the contin- 
uous blue, dashed green, dot-dashed orange and double-dot 
dashed red lines, respectively. The warm phase is reported 
in the left panel. Although showing slightly different evolu- 
tions in detail, the amount of warm gas at 2 = is almost 
the same, ~ 35 per cent, for the runs with galactic winds 
and with BH feedback. This gas should be partly responsible 
for lo cal Lyman-a forest absorption (e.g.. iDanforth fc Shulll 
2008), samples the outskirts of galactic haloes and is weakly 
sensitive to the energetics and mass-loadin g parameters of 
the e nergy-driven galactic wind feedback (jPierleoni et al.l 
l200Sl ). The case with no feedback (NW) has instead a ~ 10 
per cent lower amount of warm gas at 2 = 0. In the mid- 
dle panel, we show the evolution of the WHIM phase: with 
BH feedback the amount of WHIM at 2 — reaches 50 per 
cent, about 10 per cent more than for the runs with winds 
and 15 per cent higher than for the NW case. Similarly to 
the results for the warm phase, there are small differences 
between the W and CW runs, also in the amount of WHIM 
gas and stars. A more significant difference between W and 
CW runs is found for the mass fraction in stars. The rea- 
son for this difference is in the lower efficiency of the CW 
feedback model. Indeed, hydrodynamically coupled winds 
deposit kinetic energy through hydrodynamical processes in 
relatively higher-density environments, which are character- 
ized by shorter cooling times, which causes the thermalized 
energy to be promptly radiated away. 

The evolution of the WHIM mass fraction in the BH 
run is initially identical to that for the NW run at 2^ 3.5. 
Indeed, at high redshift there is still a limited number of DM 
haloes whose mass is large enough to host actively accreting 
BHs. At lower redshift, when accretion onto BHs is more ef- 
fective, the released feedback energy efficiently heats the gas 
surrounding galaxy-size haloes, thus moving gas from the 
warm to the WHIM phase. This explains at the same time 
the lower amount of warm gas in the BH run at intermediate 
redshift and the correspondingly larger amount of WHIM. 
As for the effect of galactic winds, they play a role in heat- 
ing circum-galactic gas already at high redshift, soon after 
the onset of star formation. This explains the larger amount 
of WHIM at z> 2.5 in the W and CW runs with respect 
to the BH run. The situation is reversed at lower redshift, 
when the depths of the forming potential wells become large, 
making gas heating with winds less effective with respect to 
BH feedback. 

In the right panel we plot the evolution of the mass 



fractions in stars and hot gas, which are represented by the 
thin and thick lines, respectively (note the different scale 
in the y— axis). As expected, the mass in stars in the runs 
with and without feedback is significantly different, largely 
reflecting the behaviour of the SFR histories shown in Fig. 
[T] the BH run has a fraction of stars which is five (two) 
times less than the NW (CW) run, as a consequence of the 
quenching of star formation taking place around 2 = 3, while 
the W run produces a slightly smaller amount of stars when 
compared to the CW one. 

Recently, iLi fc White] (|2009l ) used SDSS data to recon- 
struct the galaxy stella r mass fun c tion i n the local Universe. 
Assumi ng an IMF by IChabrierl ((2003.) (very close to the 
IMF bv lKroupall2001 I used in the simulations shown in Fig- 
ure [3), they concluded that 3.5 per cent of the baryonic 
mass is locked into stars. While clearly ruling out a model 
with no feedback, this observational estimate of the stellar 
mass fraction is close to the predictions of the BH run, al- 
though the star formation rate in this case is too high at 
3;$ z^z 6 and slightly too low at 2^ 0.5. As for the W and 
CW runs, they produce about two times and 2.5 times more 
stars than the BH run, irrespectively, as a consequence of 
a star formation excess below 2 ~ 1. These results suggest 
that both feedback mechanisms should be active: SN driven 
galactic ejecta should regulate star formation already at high 
redshift within relatively small galaxies, while BHs should 
quench cooling in large haloes around the maximum of the 
star formation, thereby keeping gas pressurised within such 
haloes down to z = 0. 

The amount of hot gas at temperatures above 10 7 K is 
instead very similar for all the simulations and around 3 per 
cent. This is a consequence of the fact that the bulk of hot 
gas lies within groups and clusters of galaxies, whose gas 
content is marginally affected by the details of the feedback. 
We note, however, that while the amount of hot gas in the 
BH run is larger than in the other runs at high redshift, it 
drops below them at z = 1, with an overall trend that is 
opposite to that of the WHIM component. Indeed, at high 
redshift BH feedback is quite efficient in heating gas to high 
temperature, thus providing a sort of diffuse pre-heating. 
At later time, this pre-heated gas has a harder time to fall 
into group-size potential wells. To quantify this effect, we 
note that in the BH run about ~ 10 per cent of the hot gas 
lies below <5b ~ 1, while only J$ 0.1 per cent of hot gas lies 
at such low densities in the other runs. 

It is worth pointing out that although at 2 = the total 
amount of hot gas is about the same, ~ 3 per cent, for all 
feedback models, its distribution as a function of density has 
a distinct pattern in the BH run. In fact, in the runs with 
no BH feedback gas is heated mostly by the process of grav- 
itational virialization within the potential wells of galaxy 
groups and clusters. As a consequence, in those runs only 
about the 5 per cent of hot gas lies below virial overdensi- 
ties, 5b ~ 50. On the other hand, a sizeable amount of gas 
displaced by BH feedback lies outside virialised haloes, with 
~ 37 per cent of it found at 6blS 50. Furthermore, the heat- 
ing efficiency of BH feedback prevents hot gas from reaching 
densities as high as in the other runs in the central regions 
of galaxy groups and clusters (see also Siiacki et al.l 120071 : 
Bhattach arva et al.l 1200 si: iFabian et alj|2009h . As a result, 
the highest density reached by hot gas in the BH runs is 



© 0000 RAS, MNRAS 000, 000-000 



8 Tornatore, Borgani, Viel & Springel 




Figure 3. The fraction of mass in the warm (left panel), WHIM (middle panel) and stellar plus hot phases (right panel, thin lines for 
stars and thick lines for hot gas) as a function of rcdshift for the NW37,256> W37,256i CW37,256 and 131137,256 runs (continuous blue, 
dashed green, dot-dashed orange and double-dot dashed red curves, respectively). 




redshift redshift 



Figure 4. Effects of resolution and box-size on the amount of gas in different phases. The fraction of mass in the warm (left panel), 
WHIM (middle panel) and stars plus hot components (right panel, thin lines for stars and thick lines for hot gas) as a function of redshift 
for the W37 i 256i W37 ! 4oo ! W75 i 5i2 runs (continuous blue, dashed green, dot-dashed orange curves, respectively). 



lower by about a factor three than in the NW, W and CW 
runs. 

Interestingly, hot metals follow the fate of the gas they 
are associated with. Indeed, the fractions of hot metal mass 
lying at different densities are quite similar to the corre- 
sponding fractions of hot gas mass. 

In order to verify the robustness of our results against 
box size and numerical resolution, we compare in Figure [4] 
the results of the W37,256 run with those of the W75.512 and 
W37,4oo runs. In general, we find that there is a very good 
convergence in the values of the mass fractions associated 
to the WHIM, warm and star phases, at all considered red- 
shifts. We only note that the amount of gas in the hot phase 
significantly increases in the W75,si2 run at zS; 1.5. This is 
due to the fact that a larger box can accommodate a larger 
number of relatively more massive galaxy groups and clus- 
ters where a larger amount of gas is shock-heated to a tem- 
perature above 10 7 K. A larger box, with size ~ 200 /i _1 Mpc, 
would certainly allow to better sample the high end of the 
halo mass function and, therefore, to provide a fully con- 
verged estimate of the mass fraction in hot gas. We note, 



however, that our larger simulation box already allow us to 
sample the scale of galaxy groups, which contain gas at a 
temperature ^ 10 7 K. Indeed, the largest halo found within 
the 75 /i _1 Mpc box has a mass of about 3 x 10 14 /i _1 Mo and 
a X-ray emission-weighted temperature of about 4 x 10 7 K. 
Therefore, we do not expect our baryon inventory in sim- 
ulations to be significantly affected by effects of finite box 
size. In the following, we will restrict our analysis to the four 
(37, 256) runs, based on the four feedback schemes (NW, W, 
CW, BH). For this reason, and unless otherwise specified, we 
omit from now on the box-size and number of the particles 
when referring to the analysed simulations. 

In summary, we conclude that at z = the amount of 
mass in the WHIM in the simulation including BH feedback 
is about 10 per cent larger than in the simulations including 
galactic winds, and also displays a different redshift evolu- 
tion. Quite remarkably, thi s result is in quantitat ive agree- 
ment with that reported bv lCen fc Ostr ikcr (2006). The fact 
that comparable values for the WHIM mass fraction are 
found with different simulation codes, based on different hy- 
drodynamic schemes, and using different implementations of 
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Figure 5. Ratio between the average age of WHIM enrich- 
ment and the age of the universe as a function of redshift 
(*WHIM Acosmic)- Results are shown for the W 37 ,256, NW 37i2 56, 
CW37 i 256 and BH37 256 simulations, represented by the continu- 
ous blue, dot-dashed green, double-dot-dashed orange and dashed 
red lines, respectively. 



feedback processes highlights that this should be considered 
as a robust prediction of models of cosmic structure for- 
mation. Therefore, if future observations will falsify these 
predictions, this will have direct implications on the need to 
include new physical processes in simulations. 



3.3 The history of enrichment 

In Figure [S] we show the redshift dependence of the aver- 
age age of enrichment of the gas particles in units of the 
age of the universe at the same redshift. Results are shown 
for warm and WHIM phases. Besides showing results for all 
gas particles belonging to each of these two phases, in the 
following we will also classify gas within each phase accord- 
ing to its overdensity. Specifically, we will define a collapsed 
phase, which is made up by all gas particles having over- 
density <5b > 50 (we note that this is the typical overdensity 
reached within virialised haloes), and a diffuse phase, which 
contains gas at Sb < 50. 

Furthermore, we define the average age of enrichment 
of a gas particle at redshift z as 



i{z) 



J2i Amz,i(z)U 
mz,i{z) 



(2) 



where the sum is over all the time-steps performed until red- 
shift z, Amz.i is the mass in metals received by the particle 
at the i-th simulation time-step, ti is the cosmic time at that 
time-step and mz{z) is the total metal mass received by a 
particle before z. According to this definition, large values of 
enrichment age, at a given redshift, correspond to a smaller 
look-back time for the epoch of enrichment, while smaller 



values of i(z) indicate a more pristine metal enrichment. In 
the limit in which all the metals were received by the par- 
ticles at the considered redshift z, then the enrichment age 
coincides with the cosmic age at that redshift. The mean of 
the ages of enrichment of a given phase is then computed by 
averaging eq. © over all particles belonging to that phase, 
each particle being weighted according to its mass in metals: 



phase \ 



J2jtj(z)mz,j(z) 



(3) 



It is worth pointing out that this way of computing the en- 
richment age of a gas phase does not provide a measurement 
of the average time of enrichment of the phase to which the 
particles belong; instead, it is an estimate of the average 
time at which the metals that belong to that phase were de- 
posited in the gas particles currently making up that phase. 
Figure [5] shows the redshift dependence of the age of en- 
richment: the left panel shows the enrichment epoch of the 
WHIM phase in units of cosmic time at the same redshift, 
i p hase{z) /tcosmic (z) , for the same simulations of Figure [2] 

Figure [5] shows the enrichment epoch of the WHIM 
phase. We note that it decreases with time in all runs, thus 
implying that at high redshifts the metal mass deposition 
takes place at a larger fraction of the cosmic time. Thus, 
the lower the value of t p hase(z) Acosmic (z), the earlier the 
enrichment of the particles belonging to that phase hap- 
pened. We also point out that the rate at which the ratio 
frwHiM(.z)Acosmic(z) decreases with redshift provides a mea- 
sure of when the bulk of metals are received by the parti- 
cles belonging to each phase. For the NW simulation, the 
WHIM enrichment is always quite recent and is due to the 
combined action of star formation and structure formation: 
while the former sets the time-scale at which metals are 
produced and assigned to the gas surrounding star-forming 
regions, the latter determines the epoch at which potential 
wells are sufficiently developed for gas dynamical processes 
to become effective in removing gas from such regions dur- 
ing the hierarchical assembly of structures. The W and CW 
runs exhibit a redshift dependence that is quite similar to 
that of the NW run. This is not surprising since the wind 
mass-load that carries the metals to the WHIM phase is also 
proportional to the star formation rate, while gas-dynamical 
processes dominate at low redshifts, when the star formation 
rate drops at a similar rate in all these runs. A significantly 
different history of enrichment is instead found in the BH 
run. In this case, the sudden episodes of energy release, oc- 
curring around 2^ zj; 4 when the mass accretion rate onto 
BHs peaks, are very effective in heating and expelling a large 
amount of enriched gas from the galaxies. Indeed, the faster 
decrease of the iwHiM(z) Acosmic (z) ratio for the simulation 
with BH feedback is a consequence of the different timing of 
the enrichment episodes of the particles that then ended in 
the WHIM. Therefore, while winds produce an enrichment 
which starts earlier and proceeds more gradually, AGN-BH 
feedback heats and enriches the particles at later epochs and 
over a shorter time interval. On the other hand, the lack of 
feedback capable to displace enriched gas in the NW run 
delays the release of metals from star-forming regions and 
closely links the timing of enrichment for WHIM and warm 
phases to the timing of growth of large potential wells, where 
stripping of enriched gas from merging haloes takes place. 
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Figure 6. Left panel: difference, in Gyr, between the average ages of enrichment for WHIM and warm particles in the diffuse phase 

(*WHIM,<5 b <so -*warm,<5 b <5o)' Ri § ht P anel: the same as in the left P anel but for § as in the collapsed phase (f WH i M ,<5 b >so -*warm,<5 b >5o)- 
Results are shown for the W37 i 256 J NW37,256i CW37 i 256 and BH37256 simulations, represented by the continuous blue, dot-dashed green, 
double-dot-dashed orange and dashed red lines, respectively. 



Another useful diagnostic to study the enrichment tim- 
ing of the IGM is represented by the relative delays of en- 
richment between the WHIM and warm phases and between 
the diffuse and collapsed environments. We recall that we 
defined a gas particle to belong to the diffuse phase when- 
ever it has density contrast 5 b < 50, while it belongs to the 
collapsed phase for 5 b > 50. 

In the left panel of Figure |5] we plot i wmM <5 b <so — 
i a , tn , while the right panel shows E,™,, s ^ rn — 

warm,O b <50' to " WHIM,0 b >50 

i j . „. As shown in this panel, at high redshifts the 

warm,O b >50 ^ ' ° 

enrichment of the diffuse WHIM phase in the NW run is on 
average roughly coeval to that of diffuse warm phase, while 
it becomes progressively more recent at lower redshifts (by 
~ 0.2 Gyr at 2 ~ 2 and by ~ 1 Gyr at 2 = 0). As already 
mentioned, this evolution is purely driven by gas-dynamical 
processes which bring metals from merging haloes to the dif- 
fuse medium, whenever these haloes enter in the pressurized 
medium permeating large potential wells. Since enriched gas 
is extracted from galaxies and groups of galaxies at temper- 
atures typically hotter than that of the IGM, cooler diffuse 
enriched particles need more time to reach lower tempera- 
tures. This implies that they have been stripped earlier than 
hotter gas, thus explaining why warm diffuse gas has been 
enriched before WHIM diffuse gas. 

Quite interestingly, the presence of galactic winds (W 
and CW runs) is effective in improving the mixing between 
the warm and WHIM diffuse phases, thus making their en- 
richment coeval down to lower redshifts, z ~ 1. At lower red- 
shifts, gas stripping within large structures becomes again 
the dominant process, thus recovering the same qualitative 
behaviour as in the NW run. As for the run with BH feed- 
back, its behaviour is quite close to that of the NW run at 



high redshifts. At 2^$ 3 gas, displacement from already en- 
riched haloes to the WHIM phase becomes gradually more 
efficient, thus making the enrichment of the diffuse WHIM 
more recent than for the diffuse warm gas. 

As for the gas in the collapsed phase (<5 b > 50, right 
panel) , the enrichment of the warm and of the WHIM phases 
is almost coeval down to 2 ~ 2 in the NW run. The reason 
for this lies in the relatively short time scale over which warm 
particles within haloes are heated to T > 10 5 K and enriched 
when they approach star forming regions. At lower redshift, 
the growth of progressively larger structures makes gas in a 
shock-heated phase to have a progressively higher temper- 
ature, until it eventually reaches the WHIM temperature 
range. At z = this gas makes up a phase with T > 10 5 K, 
reaching overdensity of up to <5 b ~ 10 4 , as shown in the up- 
per left panel of Fig. [5] Differently from the warm dense gas, 
the shock-heated WHIM does not lie close to star forming 
regions and, as such, it did not experience recent enrich- 
ment episodes. By 2 = 0, the enrichment age of the WHIM 
takes place on average ^ 1 Gyr earlier than that of warm 
gas. As for the W run, we note that the effect of winds is 
that of shortening the time scale that warm gas takes to 
reach WHIM temperatures, thereby making even more co- 
eval the enrichment age of the diffuse and collapsed phases. 
The run with BH feedback is again very close the NW run 
at z^j 3. At lower redshift, BH accretion becomes efficient. 
The subsequent energy feedback causes a fast removal of re- 
cently enriched gas from the very dense warm phase, with 
5b > 10 4 , surrounding star-forming regions. This effect is 
visible in Fig. comparing upper left and lower right pan- 
els one clearly notices the depletion of dense and warm gas 
in the BH run. This gas is shock heated to larger tempera- 
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tures, thus providing a supply of recently enriched medium 
to the WHIM phase. 

In general, these results show that the timing of en- 
richment of the WHIM does depend on the nature of the 
feedback included in the simulations. The presence of winds 
leaves its fingerprint in the timing of enrichment of the dif- 
fuse component, with <5b < 50, at relatively high redshift, 
z^j 1. On the other hand, BH feedback has a much more ev- 
ident effect at relatively low redshift, z'S 2, i.e. in the regime 
where it reaches the peak of efficiency in quenching star 
formation and in displacing gas from star forming regions. 
Quite interestingly, the effect is more pronounced for the 
collapsed gas phase, <5b > 50. This suggests that studying 
the enrichment pattern of galaxy clusters and groups, out 
to their outskirts, should represent the best diagnostic for 
the role played by A GN feedback in det ermining the cosmic 
cycle of metals (e.g.. iFabian et al1l2009t ). 



3.4 Properties of the WHIM in the local Universe 

3.4-1 Total and metal mass distribution as a function of 
density and temperature 

While in the previous sections we analysed the evolution of 
different phases, defined according to temperature criteria, 
we present now an analysis of the metal and gas content 
of the WHIM at z — as a function of gas density and 
temperature. As such, this analysis has implications for the 
detectability of the WHIM in the local Universe and for 
connecting its properties to the feedback mechanisms which 
determine its density and metallicity structure. The results 
are presented in Figures [7J and [8] for the density and temper- 
ature distributions, respectively. In the left panels of these 
figures we show the probability distribution function (PDF) 
for the total WHIM mass, while in the right panels we show 
the corresponding PDF for the metal mass in the WHIM. 
These PDFs are defined so that they provide the fraction of 
WHIM gas mass (metal mass) contributed by unit logarith- 
mic intervals in gas density (Fig. [Jj and temperature (Fig. 

Consistently with the phase diagrams shown in Fig. [2] 
the left panel of Fig. [Jj shows that BH feedback predicts a 
mass distribution for the WHIM component which is quite 
different from the other simulations: it reaches smaller un- 
derdensities, while the high-density tail is suppressed. The 
reason for this is that this feedback model provides episodes 
of strong heating of the gas around BHs, allowing it to reach 
low density regions, thus reducing the gas content of viri- 
alised haloes. At z = 0, we find that ~ 50 per cent of the 
WHIM in the NW, W and CW runs has 5 h < 50, thus ly- 
ing outside collapsed structures. This fraction increases to 
^ 80 per cent in the run with BH feedback. Correspond- 
ingly, the peak of the WHIM density distribution in the BH 
case lies at 8b — 10, a value that increases to 8b — 80 in the 
other simulations. This re sult i s in line with that fou nd by 
iBhattacharva et all l|200Sfl and lPuchwein et alj (|2008Tj . who 
showed that BH feedback is effective in removing gas from 
virialized haloes having masses typical of galaxy groups. 
Differences between the W, CW and NW runs are instead 
more prominent in the high density tail of the distribution, 
8b > 10 4 , where they involve in any case only a tiny frac- 
tion of the total WHIM mass. In the NW and CW runs, gas 



accumulates around star forming regions, reaching densities 
close to the threshold density for the onset of star formation. 
The effect of the decoupled winds in the W run is that of 
efficiently removing this gas and bringing it to lower density 
regions, with 5 h ~ 10 3 - 10 4 . 

Similar features that describe the WHIM mass distri- 
bution are also visible in the distribution of the mass in 
metals (right panel of Fig. [JJ. Metals are more likely found 
at overdensity Sb ~ 10 in the BH simulation, while for the 
NW case this overdensity is usually between 10 2 and 10 3 . 
Again, this comparison clearly shows that BH feedback is 
highly efficient in displacing metal-enriched gas from high 
density regions, thereby providing a more uniform enrich- 
ment of the WHIM at z = 0. The PDFs for the CW and 
W runs have a sort of double-peaked shape: one peak is 
at moderate overdensity, Sb ~ 10, that corresponds to the 
WHIM lying outside virialized haloes, where winds typically 
deposit metal-enriched gas; the second peak is at 8b ~ 10 3 , 
well inside virialized haloes and close to star-forming regions. 
Such a double-peak structure is more evident for the W run, 
where winds are more efficient in escaping high-density re- 
gions. Finally, we note that in the BH case the metal mass 
fraction at 8b > 10 2 is dramatically smaller (by up to 2 
orders of magnitude) than in all the other cases. 

As for the temperature distribution of the WHIM, 
shown in the left panel of Fig. [8] we find that winds have a 
negligible effect. For BH feedback, it has the effect of reduc- 
ing the mass fraction of gas above 10 6 K, while correspond- 
ingly increasing the fraction of gas at lower temperature. 
Indeed, BH feedback is efficient in removing gas from haloes 
having virial temperature within the WHIM range. After 
being heated by BH feedback, this gas leaves the halo po- 
tential wells and cools to lower temperature, T < 10 K, by 
adiabatic expansion. 

In general, the left panels of Figs.[7jand[8]highlight that 
the presence of galactic winds has a negligible impact on the 
distribution of gas mass in both density and temperature. 
Although the presence of BHs has a more apparent effect on 
the distributions of gas mass, it is clear that the distribu- 
tions of metal mass fractions are much more sensitive to the 
adopted feedback model (see right panels of Figs. [7J and |SJ). 
This is not surprising, since feedback mostly impacts on the 
gas surrounding star forming regions, which has the higher 
enrichment level. In the NW run, the metal mass fraction in 
gas at the low-temperature boundary of the WHIM phase, 

10 5 K, is about one order of magnitude lower than for 
the hot WHIM at T< 10 7 K. The effect of winds is that of 
transporting metal rich gas from high density regions, where 
it reaches high temperature, to lower density regions, where 
is cools down by adiabatic expansion. As expected, hydro- 
dynamically coupled winds (CW) make it harder for the gas 
to leave the densest regions, thus producing a distribution 
which is intermediate between those of the NW and W runs. 

3-4-2 The distribution of the WHIM metallicity 

Having characterized how the mass in gas and in metals is 
distributed as a function of the WHIM density and tem- 
perature, we describe now how the WHIM metallicity de- 
pends on density and temperature. This characterization of 
the WHIM is more observationally oriented. Indeed, differ- 
ent distributions of metallicity determine what are the best 
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Figure 7. Probability distribution functions (PDF) of the total WHIM mass (left panel) and WHIM mass in metals (right panel) at 
z = as a function of gas density (in units of the mean cosmic baryon density) for the W, NW, CW, BH simulations, represented by 
the continuous blue, dot-dashed green, double-dot-dashed orange and dashed red lines, respectively. 




Figure 8. Probability distribution functions of the rotal WHIM mass (left panel) and WHIM mass in metals (right panel) at z = as 
a function of gas temperature for the W, NW, CW, BH simulations, represented by the continuous blue, dot-dashed green, double-dot- 
dashed orange and dashed red lines, respectively. 



tracers (chemical elements and their ionization states) to re- 
veal the presence of the WHIM and characterize its physical 
properties: the presence of different ionization species, that 
are thought to be the best WHIM tracers, depends in fact 
upon loc al conditions of dens ity, temperature and ionization 
field (see lRichter et al1 l2008. for a review). 

We show in Figure|9]the distribution of the total metal- 
licity as a function of density within the WHIM tempera- 
ture range (left panel) and as a function of temperature 
(right panel). Since the metallicity distributions have always 



a quite large scatter, a meaningful way of presenting the re- 
sults is in terms of the mean (continuous lines), median 
(dashed lines) and of the 10 and 90 percentiles (shaded 
area for the W run, dot-dot-dashed lines for the BH run) of 
these distributions. For reasons of clarity, we do not report 
in these figures the scatter for the NW and CW run. 

As for the density dependence of metallicity in the W, 
CW and NW simulations (upper left panel), we note that 
at overdensities between 10 and 100, which is the range 
where the PDF of the WHIM mass distribution reaches 
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its maximum value (see Fig. [7}, the median metallicity has 
a strong positive correlation with gas density and increase 
from log(Z/Zo) ~ -6 to Iog(Z/Z s ) ~ -2. Furthermore, 
the overall scatter around these median values is about two 
orders of magnitude, so that metallicity values as large as 
O.IZq are not unlikely. An obvious observational implication 
of this large scatter in the metallicity distributions is that 
a fairly large number of lines-of-sight are required, along 
which measuring WHIM metallicity, in order to properly 
populate such a scatter. It is worth mentioning that the 
large difference between the average and median metallic- 
ities in regions with 8b^ 1 arises from the highly skewed 
distribution of metallicity in this regime: most of the under- 
dense gas is not enriched at all, while metals mainly lie in 
a small number of highly enriched particles with Z£ 0.01. 
This is mostly true for BH and NW runs, while in W and 
CW runs the distribution of metals with particle metallicity 
is significantly shallower. We point out that highly enriched 
particles preserve their metal content due to the intrinsic 
lack of diffusivity of the SPH . Including an ex plicit descrip- 
tion of metal diffusion fe.g.. [Greif et a.1.1 12009T ) would make 
such particles sharing their metal content with the surround- 
ing metal-poor particles. 

As for the metallicity distribution as a function of tem- 
perature (bottom panels), we note a clear trend for metallic- 
ity to increase with temperature in all cases, although the de- 
tails of this trend have a dependence on the feedback model. 
The median metallicities in the W, CW and NW runs (bot- 
tom left panel) increase from 10 -4 — 10~ 3 Zq to ^ 10 _2 Zq 
as the temperature increases from 10 6 K to 10 7 K. The NW 
simulation predicts a metallicity of the hottest WHIM, with 
T ~ 10 7 K, which is higher by about a factor three than for 
the W and CW runs. This gas is located within the virialised 
regions of galaxy group haloes, which have typical overden- 
sities of about 10 2 — 10 3 (see upper left panel). In the NW 
run, this gas, which remains at these overdensities for a rela- 
tively long time before being stripped, has been more heav- 
ily enriched by the more (in fact, exceedingly) intense star 
formation. Quite interestingly, the CW simulation shows in- 
stead larger median metallicity at temperatures 10 5 ' 5 — 10 6 
K, with a shallower decline towards lower temperatures. Gas 
at this temperature is a mixture of highly enriched gas in 
virialized haloes and of poorly enriched gas at low density. 
As can be inferred from the phase diagrams shown in Fig. 
OH gas at z = in the CW run lies within the WHIM tem- 
perature range either for 10$; 8b $ 100 or right outside of 
the star-forming regions, i.e. 10 4 ^ <5b ^ 10 5 , where enriched 
wind particles loose their kinetic energy due to hydrody- 
namical interactions. As a consequence, in the CW run the 
WHIM at these overdensities is more enriched than in other 
runs. 

As for the comparison between NW and BH runs (right 
panels), we note that BH feedback removes a lot of metal 
enriched gas from the surroundings of star forming regions 
(see also Fig. [2}. This explains the lower metallicity of the 
densest WHIM in the BH simulation (upper left panel). Fur- 
thermore, the suppression of star formation in the BH run 
also reduces the enrichment level of the shocked gas within 
collapsed regions with 10 2 < 8 h < 10 3 . At the same time, BH 
feedback is quite effective in increasing the metallicity of gas 
lying outside virialized structures, down to rather low den- 
sities, 8 b ^ 10. The reason for this widespread enrichment 



is the high efficiency of BH feedback to displace highly en- 
riched gas from the haloes of large galaxies in the redshift 
range z ~ 2-4, in correspondence of the peak of the BH 
accretion rate. 

As for the dependence of the WHIM metallicity on tem- 
perature, the NW and BH runs provides similar trends, with 
a lower level on enrichment in the BH run. 

In general, the results of our simulations confirm that 
different feedback models leave distinct signatures in the 
metallicity structure of the WHIM. While these results 
have in principle interesting implications for the observa- 
tional characterization of WHIM, it must be emphasized 
that the large scatter expected in the metallicity distribu- 
tions requires that metallicity measurements should be car- 
ried out for a fairly large number of lines-of-sight. We defer 
to a future paper the presentation of an observationally- 
oriented analysis of our simulations, where we discuss the 
efficiency with which the metallicity structure of the WHIM 
can be recovered from mock absorption spectra of back- 
ground sources. 

3.5 Evolution of the relative metal abundances 

While the above analysis was aimed at describing the overall 
metallicity of the different gas phases, a much richer amount 
of information is provided by our chemo-dynamical simula- 
tions. Indeed, the possibility of tracing the enrichment of 
different elements allows us to study differences in the tim- 
ing of enrichment of different heavy elements. In this section 
we explore the redshift evolution of the ratios between the 
abundances of Oxygen with respect to Iron and to Carbon, 
with the purpose of quantifying the signatures that IMF and 
feedback mechanisms leave on the WHIM enrichment his- 
tory. Atomic transitions associated to these three elements 
are often used to reveal the presence of the diffuse IGM and 
to measure its enrichment level at different redshifts. For in- 
stance, CV, OVI and OVII lines are considered as the most 
prominent absorption features in the far-UV and soft X-ray 
spectra of background sources to re veal the presence of the 
WHIM in the nearby Universe (e.g.. iRichter et alJliooj . for 
a review). At the same time, CIII, CIV and OVI absorp- 
tion lines are common ly used to trace the metal content of 
the I GM at z^ 2 (e.g.. lSchave et"aL I l2003l ; iPieri fc Haehneld 
I2004T ). Finally emission features associated to the OVIII line 
and to the Fe-L and Fe-K complexes are used to trace the 
meta l content of the in tra-cluster and intra-group media 
(e.g.. IWerner et al"1l2008l . for a review). 

Since different elements are produced in different pro- 
portions by different stellar populations, the evolution of 
their relative abundances is expected to depend, for a fixed 
mass-depende nt life-time funct i on, on the shape of the stel- 
lar IMF (e.g., Tornator e et al. 1 120071 ; IWiersma etHI 120091 ; 
see lBorgani et aL1l2008l for a review). In general, we expect 
that different spatial distributions characterize the enrich- 
ment pattern for different elements. As an example, prod- 
ucts of SN-II are released over a shorter time-scale than 
those arising from SN-Ia. Since star particles are expected 
to move from their original location where they have formed, 
e.g. due to merging or stripping processes, we expect that 
SN-II products pollute gas particles lying very close to star 
forming regions, while SN-Ia products should have a rela- 
tively more diffuse distribution. We expect this effect to be 
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Figure 9. Total mctallicity of the WHIM (gas particles at temperatures 10 5 — 10 7 K) in solar units at z = 0, as a function of gas 
density (in units of the cosmic mean baryon density (pt> )> left panel) and temperature (right panel). In each panel, the grey shaded area 
encompasses the 10 and 90 percentiles of the W run, while the dot-dot-dashed lines show the same percentiles for BH the run. Thick 
coloured lines show the average metallicities while thin dashed lines show the median metallicities. 



more apparent in the dense environment of galaxy clusters 
where a population of diffuse inter- galactic stars is gener- 
ated by dynamical processes. Indeed. rTornatore et alJ (|2007l ) 
found from chemo-dynamical simulations of galaxy clusters 
that the distribution within the intra-cluster medium (ICM) 
of the metals produced by SN-II is more clumpy than that 
provided by SN-Ia, a pr ediction that has been al so confirmed 
by observational data (jSivanandam et aT]|2009l ). 

In Figure [10] we plot the [C/O] and [O/Fe] ratiofl as 
a function of redshifts for the four models and for both 
the warm (left panels) and WHIM (right panels) phases. 
These three elements are produced in different proportions 
by SN-II and SN-Ia: Oxygen is almost entirely produced 
by SN-II, and Iron is largely provided by SN-Ia; while Car- 
bon is produced in comparable proportions by SN-II and 
SN-Ia, the main contributors to it at timescale compara- 
ble to those of SN-Ia are AGB stars. The [C/O] ratio has 
been widely in vestigated in a variet y of environments, from 
dwarf galaxies l|Garnett et al.lll995), to low redshift ICM a s 
probed by QSO absorption lines ilPanforth fe Shulll200Sl) . 
to high redshift IGM at z = 2 - 4 (jAguirre et alj|2008j ). As 
for the [O /Fe] ratio, it is traced in emission with X-ray spec- 
troscopy of relatively poor g alaxy clusters and groups (e.g., 
iRasmussen fe Ponmanll2007t) and is used to study the rela- 
tive role played by SN-II and SN-Ia in enriching the ICM. 

Despite the fact that different feedback prescriptions 
induce rather different evolutions for the [C/O] ratio, espe- 
cially for the WHIM phase, the values attained at z — 
are rather similar for all the simulations. We find [C/0]~ 
-0.2 and [C/0]~ -0.15 for the warm and for the WHIM 



2 Following a standard notation, we define the relative abun- 
dance between the elements X and Y as [X/Y~] = log(Zx /Zy) — 
log(Z x ,Q/Z Y ,e)- 



phases, respectively. These values can be compared with a 
sim ilar analysis of chemo- d ynam ical simulations performed 
by lOppenheimer fe Pavel feoOSl) . At z = they found 
[C/OJ~ -0.09 and [C/OJ~ -0.05 for the warm and warm- 
hot phases, respectively, thus about 0.1 dex higher than our 
results (note that we rescaled their Anders & Grevesse 1989 
abundances to o ur used values). A possible r eason for this 
difference is that lOppenheimer fe Pavel (120081 ) used a model 
for momentum-driven galactic ejecta. However, since they 
obtain rather stable results for a variety of feedback schemes, 
it is not clear whether this is a likely explanation for the 
difference. Another possibility lies in the different sets of 
yields, mass limit for the SN-II progenitors (10M©, instead 
of 8M© as in our simulatio ns), and the mass-dependen t 
life-time function adopted by lOppenheimer fe Pavel (2008). 
For instance, assuming a larger limiting mass for SN-II 
turns into a relatively lower amount of Oxygen produced, 
which could explain t he la rger value of [C/O] found by 
lOppenheimer fe Pavel l|2008l ). 

As for the evolution of the WHIM phase in the BH run, 
we note that [C/O] is very close to the NW values at high 
redshift. This is expected, since at early times BH accretion 
is quite ineffective. At these high redshifts, gas at WHIM 
temperatures starts involving shock-heated gas in filaments. 
This diffuse gas has been relatively more enriched by long- 
lived stars that had time to move away from the highest 
density star-forming regions. This explains why the relative 
abundance of Carbon increases with respect to Oxygen, a 
trend that extends in the NW run down to z ~ 2, when a 
nearly solar relative abundance is attained. After this red- 
shift, the diffuse gas whose temperature reaches the WHIM 
values, has a progressively lower degree of enrichment. At 
the same time, the large amount of low-redshift star forma- 
tion in the absence of feedback causes a prompt release of 
Oxygen from short-lived stars, thus motivating the gentle 
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Figure 10. Evolution of the mean value of the [C/O] (upper panels) and [O/Fc] (lower panels) ratios for the warm (left panels) and 
WHIM (right panels) phases. Each panel shows the results for the W256,37i NW256,37, CW256,37 and BH256.37 simulations, represented 
by the continuous blue, dashed green, dot-dashed orange and double-dot dashed red lines, respectively. 



decline of the [C/O] ratio. As for the run with BH feed- 
back, the WHIM starts receiving at z ~ 3.5 a significant 
contribution from gas expelled from haloes where gas ac- 
cretion onto BHs reaches the peak of its efficiency (see also 
the central panel of Fig. [3]). At lower redshift the star for- 
mation is strongly quenched in the BH run, thus explaining 
the substantial flattening of [C/O] for both the warm and 
the WHIM phases. 

The evolution of the [C/O] ratio for the CW and W 
simulations is similar for the warm and the WHIM phases, 
for which it rises gently by 0.1 dex and 0.2 dex, respectively, 
from z = 4 to z — 0. Differently from the BH feedback, 
winds start affecting the level of star formation already at 
high redshift, z > 4. This explains the slower increase of 
[C/O], which extends down to low redshift, for both the 



warm and the WHIM phases. While the W and CW runs 
provide very similar results for the warm phase, the tendency 
of the decoupled winds to transport more efficiently enriched 
gas outside haloes justifies the slightly larger values of [C/O] 
in the WHIM phase, as found in the CW simulation. 

The values obtained at z = 2 — 4 can be compare d with 
those inferred for the IGM by lAguirre et all (2008) using 
the pixel-optical depth technique in high resolution QSO 
spectra. Their findings of course depend on the strength 
and shape of the assumed UV background. For a back- 
ground which is made by galaxies and quasars they find 
[C/0]= —0.7 ± 0.2 (this value refers to Anders & Grevesse 
1989 abundances and the reported errorbar corresponds to 
systematic uncertainties), which is in reasonable agreement 
with the results of the warm phase (left panel) that better 
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Figure 11. The effect of changing the IMF on the evolution of the [C/ O] ratio for the w arm (left panel) and WHIM (right panel) phases. 
The solid blue curves are for the reference run based on the IMF by iKroupal l|200lf ) (W37,256 run), the dashed gree n curves are fo r 
the top— heavy IMF by lArimoto &; Y oshii (198^) (Wayay^se run) and the dot-dashed orange curves are for the IMF by S alpeterl dl955l ~) 
^837,256 run). 



samples the relatively low densit y IGM. From the analysis 
made bv lDanforth fc Shulll (|2008l ). using absorption lines in 
the low redshift IGM, values of [C/O] in the range [—1,0] 
can be inferred, depending on which ion is used as a tracer 
of the metallicity. This result is again in broad agreement 
with the results of our simulations. 

Since Oxygen is mostly contributed by SN-II, the be- 
haviour of the [O/Fe] ratio has the opposite trend compared 
to [C/O], with an even more pronounced redshift depen- 
dence. The decreasing trend at low redshift is due to the 
fact that Fe is primarily produced by SN-Ia that have long- 
lived progenitors. The values of [O/Fe] at z — are ~ —0.15 
for the warm phase and ~ —0.3 for the WHIM. At z > 3, 
this abundance ratio approaches the solar value. Our find- 
ings are in good agreement with the [O/ Fe] ratios of gas 
partic les in the ICM of isolated clusters of iTornatore et ail 
|2007l ). in particular values of —0.2 are reached in the outer 
parts of galaxy clusters and this is close to the value of 
the WHI M [O/Fe] ratio. As for the c omparison with the 
results by lOppenheimer fc Pavel (120081 ). we note that our 
predicted values of [O/Fe] are about 0.4 dex smaller than 
their WHIM and warm phases once we consider the differ- 
ence abundances used. These differences highlight the need 
to carry out detailed comparison between chemo-dynamical 
codes to distinguish the effect of the different implementa- 
tion of the chemical evolution model from the effect of the 
different feedback models introduced and of the prescrip- 
tions to distribute metals around star forming regions. 

As a final analysis, we verified the effect of changing 
the IMF on the resulting evolution of [C/O]. Since the 
slope of the IMF determines the frequency of different SN 
types, we expect that a top-heavier IMF provides a rela- 
tively larger number of SN-II and, therefore, a correspond- 
ingly lower value of [C/O]. This is indeed confirmed by the 
results in Figure ITTl which show a lower value of C/O. A dif- 



ference of about 0. 15 dex is found for the top-heavy IMF by 
lArimoto fc Yoshiil (|l987T l. for both the warm and the WHIM 
phase , whil e no sizab l e diffe rences exist between the lKroupal 
|200ll ) and ISalpeterl l| 19551 ) IMFs. It is interesting to note, 
while using a top-heavier IMF changes the value of [C/O], 
its redshift dependence remains essentially unchanged. As 
such, the effect of changing the IMF is different from the 
effect of changing the feedback model, which instead does 
change the evolutionary pattern of the relative abundances. 

The main conclusion of this section is that, in spite of 
the dramatic differences in the feedback schemes and star 
formation histories, the C/O and O/Fe ratios conspire to 
reach comparable values at z = for all the simulations, 
with variations of about 20-25 per cent. Larger differences, 
of up to a factor ~ 2.5, are instead found at higher red- 
shift, z = 2.5. Much smaller differences in the evolution of 
the relative abundances are instead found by changing the 
IMF. This implies that an observational determination of 
the [C/O] evolution would provide information about the 
nature of the feedback mechanism responsible for regulating 
star formation and distributing metals. 

4 CONCLUSIONS 

We have presented results from an extended set of chemo- 
dynamical simulati ons which hav e been carried out with the 
GADGET-2 code (|Springelll2005l ). includin g the implemen- 
tation of chemical evolution described by ITornatore et al.l 
(2007). The analysis presented here was focussed primarily 
on the low-redshift properties of the inter-galactic medium 
(IGM) and its evolutionary properties, by considering both 
gas in a warm phase (gas particles with temperature in the 
range 10 4 -10 5 K) and in the so-called WHIM phase (gas 
particles with temperature in the range 10 5 - 10 7 K). The 
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main purpose of the analysis was to quantify the effect that 
different feedback mechanisms leave in the pattern of metal 
enrichment. Besides a simulation not including any efficient 
feedback (NW run), we performed simulations including en- 
ergy feedback resulti ng from accretion onto super-massive 
black holes (BH run; ISpringel et al]|2005l : IdT Matteo et al.l 
l2005h . and from galactic winds powered by supernova (SN) 
explosions. We considered two implementations of the latter, 
based either on locally decoupling winds from hydrodynam- 
ics, so as to allow them to leave star forming regions (W 
runs), or on keeping gas in the winds always hydrodynami- 
cally coupled with the surrounding medium (CW runs). 

The main results of our analysis can be summarized as 
follows. 

(a) The temperature-density and metallicity-density phase 
diagrams at z = are affected by the different feedback pre- 
scriptions. Including galactic winds has the effect of enrich- 
ing with metals the low-density IGM. BH feedback is even 
more effective than winds in transporting metal- enriched gas 
from the star-forming regions to the WHIM phase. BH feed- 
back also causes the presence of a non-negligible amount of 
relatively low-density, metal- enriched hot gas (10 6 — 10 7 K) 
which is not present in the simulations with galactic winds. 

(b) The fraction of baryonic mass associated to warm 
gas, with T = 10 4 -10 5 K, which should be associated to 
UV/Lyman-a absorption systems, varies from <; 90 per cent 
at high redshift, z > 3.5, to about 30 per cent in the local 
universe, with a weak dependence on the adopted feedback 
model. This is in agreement with recent observational est i- 
mates from UV spectroscopy (e.g. iDanforth &~Sh ull 2008). 
The hot phase, with T > 10 7 K, sums up to 3 per cent at 
z — 0, again almost independent of the feedback model. The 
fraction of baryons in stars is instead the most sensitive to 
feedback: it ranges from about 10 per cent in the run with 
no feedback (NW), to about 2 per cent in run with BH feed- 
back. 

(c) The WHIM phase comprises about 35 per cent of the 
baryon budget at z = for the NW run, a value that in- 
creases to about 50 per cent for the run with BH feedback. 
This confirms that the share of cosmic baryons in the dif- 
ferent phases does depend on the assumed feedback model. 
The redshift evolution of the mass fraction in the WHIM 
also differs in the different runs, with a stronger evolution 
in the run with BH feedback. 

(d) The average age of enrichment of the warm and of the 
WHIM phases differs in the different runs, by an amount 
which depends on the gas overdensity 8b- Diffuse gas (<5b< 
50) in the warm phase at z = is typically enriched 0.5-1 
Gyr earlier than in the WHIM phase. As for the gas within 
"collapsed" regions, with 5^, > 50, the warm phase is enriched 
1 Gyr after the WHIM phase in all runs not including BH 
feedback. In the BH run, enrichment of the densest WHIM 
as measured at z = takes place more than 1.5 Gyr after 
that of the warm medium. BHs enrich the WHIM more 
promptly than W in the redshift range 2< z< 4. On the 
other hand, winds are more effective at higher redshift, when 
BH accretion is still inefficient, and below z ~ 2, when star 
formation in the BH run has been quenched. At all epochs, 
the hydrodynamically coupled winds (CW run) provide a 
slightly more recent enrichment of both warm and WHIM 
phases when compared to hydro-decoupled winds (W run). 



BH feedback provides a faster enrichment at z > 2, while 
below this redshift winds provide a more prompt WHIM 
enrichment. In particular, the model with hydrodynamically 
coupled winds (CW) provides at z < 3 the most prompt 
enrichment of both warm and WHIM phases. 

(e) In order to address the multi-phase nature of the 
WHIM, we compute the distribution of gas and metals in 
the WHIM phase as a function of overdensity. As a result 
of the strong heating provided by BH feedback at z ~ 2-4, 
the characteristic density of the WHIM in the BH run is 
a factor of a few lower than in the other simulations and, 
correspondingly, a smaller amount of gas is present in the 
dense WHIM. The same trend is also visible for the amount 
of metals present in the WHIM. Typically, most of the met- 
als lie at over-densities between a few and 10 in the run with 
BH feedback, while the typical density of the metal-enriched 
gas is about two orders of magnitude higher in the run with 
galactic winds. The reason for this is that most of the metals 
are ejected by BH heating from galaxies between z — 4 and 
z — 2. The high entropy level reached by enriched gas heated 
by BH feedback makes hard for it to be re-accreted within 
collapsed haloes at lower redshift. In turn, the drop of star 
formation in the BH run causes a comparable suppression 
of metal production at z < 2. 

(f) Underdense WHIM regions (voids) have a very low me- 
dian metallicity, however the metallicity of individual gas 
particles can reach values larger than 0.1 Zq within the 90 
percentile. These enriched particles have been trasported to 
underdense regions both by galactic ejecta and by the action 
of gas-dynamical processes. Due to the intrinsic lack of dif- 
fusion in SPH, such particles spuriously preserve their high- 
metal content, rather than diffusing metals to surrounding 
metal poor particles. 

(g) The values of the [C/O] and [O/Fe] relative abundances 
at z — are similar for both the WHIM and warm phases, 
also irrespective of the feedback model. Their evolutionary 
pattern is instead sensitive to the adopted feedback model. 
Using a top-heavier IMF decreases the value of the [C/O] 
ratio by about 0.15 dex at all redshifts. 

The simulations presented in this paper have been anal- 
ysed to study separately the effects that SN-triggered winds 
and BH energy feedback have on the evolution of the IGM. 
Clearly, in a realistic situation one expects SN and AGN 
feedback to be both at work and to cooperate in determin- 
ing the cosmic cycle of baryons. However, we point out that 
our analysis was not aimed at establishing a best-fit feed- 
back model, which is able to reproduce a variety of observa- 
tional results. We aimed insted at determining the imprints 
that different feedback models leave of the properties of the 
IGM and their possible observational signatures. Further- 
more, it is worth reminding that the parameters that we 
adopted for winds and BH feedback have been fixed by re- 
quiring each of these two feedback sources to reproduce a 
specific observ ational constraint: the cosm ic star formation 
rate for winds dSpringel fc Hernquistll2003|) and t he Mbh~o 
relation for BH feedback |Pi Matteo et al.l (|2005T ). Once the 
two mechanisms are allowed to be both present in a simula- 
tion, they are expected to have non-trivial interplays, which 
necessarely require a re-calibration of their characteristic pa- 
rameters. 

The results obtained from our analysis have interesting 
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implications for the detectability of the WHIM and for the 
possibility of characterising it s thermal and chemical prop- 
erties (e.g.. lStocke et aLll2007l ). Indeed, the possibility of de- 
tecting the WHIM through absorption lines in the spectra of 
background sources does not only depend on the amount of 
mass in this phase, but also on how this mass is enriched and 
distributed in density and temperature. With our analysis 
we have demonstrated that such distributions are rather sen- 
sitive to the adopted scheme of energy feedback. Analyses 
aimed at discussing the WH IM detectability from simula- 
tions, both in e mission (e.g. 



absorption (e.g. JCen et al 



' Yoshikawa et alj 20041 ) and in 
200 it I Viel et al.ll2005h have been 



so far based on approximate descriptions of the pattern of 
chemical enrichment. We will present in a future paper an 
observationally oriented analysis of our simulations, aimed 
at quantifying how the WHIM properties can be recovered 
under realistic observational conditions in the presence of 
different feedback schemes. There is no doubt that the pos- 
sibility of performing high-resolution spectroscopy both in 
the X-ray (e.g., with micro-calorimeters onboard of large col- 
lecting area X-ray telescopes) and in the UV band (i.e., the 
now operating Cosmic Origin Spectrograph onboard of Hub- 
ble Space Telescope) will provide a leap forward in the study 
of the diffuse warm-hot baryons in the local universe. Cos- 
mological hydrodynamical simulations, like those presented 
in this paper, offer the natural interpretive framework for 
these future observations, which will not only complete the 
census of baryons at low redshift but also characterize their 
physical properties. 
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